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Abstract 

We use geodesic probes to recover the entire bulk metric in certain asymp- 
■ totically AdS spacetimes. Given a spectrum of null geodesic endpoints on the 

\ boundary, we describe two remarkably simple methods for recovering the bulk 

' information. After examining the issues which affect their application in prac- 

\ tice, we highlight a significant advantage one has over the other from a com- 

0^ ' putational point of view, and give some illustrative examples. We go on to 

\ consider spacetimes where the methods cannot be used to recover the complete 

bulk metric, and demonstrate how much information can be recovered in these 



cases. 



Ph. 1 Introduction 

The holographic principle has inspired many ways of exploring different spacetime 
geometries. The basic idea of holography, that physics in a region of space can be 
^ ' described by the fundamental degrees of freedom on its boundary, was originally 

■ applied to the area of quantum gravity by 't Hooft P and Susskind |2j, but it 

was Maldacena jSl who developed it into something more tangible: the AdS/CFT 
correspondence jlj (For a comprehensive review, see 0). 

Maldacena's conjecture postulated that string theory with AdS^ x boundary 
conditions is equivalent to = 4 super Yang-Mills theory, and AdS-type spacetimes 
have become the subject of much work over the past decade. Many authors |7| |H1 
El lin| 111! I12| I13| 1141 among others] have examined how to extract information from 
behind the horizon of black holes using the boundary correlation functions of the 
CFT. From the spacetime point of view this involves using geodesies (both null and 
spacelike) travelling from one boundary to another as probes. However, this focus 
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on singularities leaves open the question as to how much of a non-singular spacetime 
can be recovered using geodesic probes.^ 

What we aim to explore in this paper is how accurately and completely one can 
use the data contained within the boundary correlators to recover the bulk metric 
in a static, spherically symmetric, asymptotically AdS spacetime. We examine a 
variety of different scenarios, and rather than merely looking at whether a recovery 
of the metric is possible, we investigate the factors affecting its recovery in practice. 

We make use of the relation between null geodesic endpoints and singularities 
in the correlation functions of operators inserted on the boundary. In other words, 
that the correlation function between any two points on the boundary which can 
be connected by a null geodesic path diverges. Thus from the field theory we can 
identify, say, the set of endpoints of null probes which all originated from the same 
point on the boundary. This endpoint data can then be used to recover information 
about the bulk metric. 

In this paper we develop two methods for extracting this information, and 
the essence of both is that one can focus on specific radii by systematically us- 
ing geodesies which penetrate to different depths in the bulk. For null geodesies, 
there is only one effective parameter which determines the minimum radius obtained 
by the geodesic in a given spacetime, namely the ratio of the angular momentum to 
the energy. When this ratio is close to one, the geodesic probe travels around the 
edge of the spacetime; as it is reduced, the probe travels deeper into the spacetime 
before returning out to the boundary. 

What we discover is that in certain types of metric, the process of using null 
geodesic probes which penetrate to different depths can be used to recover the 
entire bulk metric in a remarkably simple manner. 

We begin in section 2 by briefly reviewing some background material on geodesies 
and their behaviour in a five dimensional AdS spacetime. We observe that in this 
case, null geodesies always travel between antipodal points on the boundary, as 
shown in Figl^j 

In section 3 we move on to considering modified spacetimes which are still asymp- 
totically AdS. In these metrics, null geodesies no longer necessarily terminate at 
their antipodal points after travelling through the bulk; we are presented with a 
spectrum of endpoints, which provides the data we need to recover the metric. We 
discover that on a plot of the endpoints (see FigElfor an example), the gradient at 
any point is related to the ratio of the angular momentum and the energy of the 
corresponding null geodesic. Using this relation together with the endpoint data, 
we devise a simple iterative method for extracting the bulk metric, beginning close 
to the boundary where the spacetime can be taken to be approximately AdS. After 
considering whether this idea can also be applied in a more general scenario, we go 
on to develop a second method for iteratively recovering the metric. 

These two methods are analysed in more detail in section 4, where we examine 

^For more general work on extracting behind the horizon information and holographic represen- 
tations of objects in AdS spacetimes see e.g. [T^ITHIITTIITSIIT^I^ETII^ . 
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a number of issues arising from using them in practice. We propose several possible 
resolutions, in each case attempting to strike a balance between the accuracy of 
the recovered metric information and the computational effort required to do so. 
We then describe a modification to the second method which greatly improves the 
accuracy without a corresponding increase in computational effort. Examples of the 
methods being used to recover the metric information in two different spacetimes 
are then given in section 5, where the superiority of the modified second method 
becomes apparent (see figures IHHHl and the corresponding tables). 

Section 6 looks at limitations in the methods which arise not from the numerics 
but from the geometry of certain spacetimes we would like to consider. We examine 
how much of the bulk metric can be recovered in setups where the metric gives a 
non-monotonic effective potential for the geodesies, and where there is a singularity 
at the centre of the spacetime. 

In section 7 we discuss the results and mention possible extensions of the project, 
such as using the endpoint data of spacelike geodesies or considering non-static, non- 
spherically symmetric spacetimes. 



2 Background 

The metric for a five-dimensional Anti-de Sitter spacetime is often given in the form: 

ds' = -f{r)dt' + ^+r'dnl (1) 
fir) 

fir) = 1 + ^ (2) 

where R is the AdS radius. 

To calculate the geodesic paths, we first suppress two of the angular coordinates 
by setting the geodesies to lie in their equatorial plane. Then from the form of the 
metric we can see the existence of two Killing vectors, d/dt and d/d(j), where (p is 
the remaining angular coordinate, which lead to two constraints on the motion. The 
first constraint is related to the energy and the second to the angular momentum: 

E = f{r)i (3) 



L = r^(j) (4) 



where = for some affine parameter A. 

The geodesic paths can be found by extremizing the action 



S= / ±'d\ (5) 
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Figure 1: Plot of the effective potential in pure AdS (with R = 1) for a sample of 
nuh geodesies, with E = 1.0 and L = 0.0,0.1,0.2,0.4,0.6 and 0.8 

which leads to a third constraint on the motion, given by: 

x^ = -f{r)f + ^^+r^^ (6) 

where x"^ = k = +1, —1,0 for spacelike, timelike and null geodesies respectively. 

By substituting Q and © into these three equations can be combined 
together to eliminate t and (j), giving: 

which can be rewritten so as to introduce the concept of an effective potential for 
the geodesies: 

r^ + yeff = ^ (8) 

with 

K// = -(/WK + i?^-^) (9) 

Note that the RHS of is zero, and so any part of the effective potential which 
is positive represents a potential barrier, see FigQ The minimum value of r obtained 
by spacelike and null geodesies incoming from r = oo is given by the largest solution 
to Veff = 0, and is the endpoint (starting point) of the ingoing (outgoing) part 
of the geodesic. The paths of geodesies can be calculated and plotted numerically 
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using these equations, and for null geodesies in pure AdS, they can also be calculated 
analytically (see appendix A). 

For null geodesies in AdS space, their endpoints always lie at the antipodal point 
no matter what the ratio of angular momentum to the energy. They can be plotted 
using the compactification given in appendix A, which rescales the boundary at 
r = oo to lie on a circle of radius 7r/2, see Fig|2l 

As mentioned in the introduction, the endpoints of null geodesies correspond 
to singularities in the set of correlation functions in the field theory lying on the 
boundary of the spacetime. For the rest of the paper we assume that the endpoint 
data we use can be generated in such a fashion from the field theory. 

3 Asymptotically AdS spacetimes 

We can consider a small modification to the pure AdS spacetime described by 
and (j2I), by replacing with: 

2 

f{r) = l + ^-pir) (10) 

where p{r) is an analytic function which behaves as for small r and r^^ for large 
r, such that the metric is non-singular everywhere.^ 

3.1 Analysing the geodesic endpoints 

In this modified spacetime, null geodesic paths do not always end at their antipodal 
point on the boundary; their final t and (j) coordinates depend on the ratio of the 
angular momentum (L) to the energy (E). From a physical point of view this can 
be thought of as the extra term p{r) representing an attractive modification to the 
centre of the spacetime, such that the geodesies follow a more curved path through 
the bulk, and thus "overshoot" the antipodal point. Their extra path length then 
also accounts for the increase in the time taken to reach the boundary at infinity. 

If we consider the set of null geodesies all beginning at the same point, (to, (po) 
on the boundary, we can obtain the full spread of endpoints by varying L/E from 
to 1. FigEl shows an example of the geodesic paths in a modified AdS spacetime 
alongside a plot of their endpoints. 

From this plot of the endpoints, what information about the bulk metric can 
be recovered? For a geodesic beginning at (to,0o) and ending at (ti,0i) on the 
boundary, we have: 

/ dt = 2 dr (11) 

J to Jrmin f{r)Jl - y'^^ 

^The behaviour of p{r) at small r ensures that we avoid a conical singularity at the origin. 
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Figure 2: A plot of some null geodesic paths passing through an AdS spacetime, all 
starting from the same point t = 0, = on the boundary at r = cxd (compactified 
to lie at arctan(oo) = ^ in the diagram). 



6 



Figure 3: A plot of the null geodesic paths passing through a modified AdS space- 
time, all starting from the arbitrary point i = 0, = on the boundary and 
with positive angular momentum. The corresponding full spectrum of null geodesic 
endpoints for this spacetime is shown on the right. 



dr (12) 



where r^m is minimum radius obtained by the geodesic, and we have introduced 
the parameter y = L/E to denote the ratio of the angular momentum to the energy. 

We note that for a fixed metric, the final time and angular coordinates will be 
functions of y only, and so we can write: 

ii 

dt = h-to = t{y) (13) 

to 

d(/. = 01 - 00 = Hy) (14) 

<f>0 

If we define the function g{y,r) as: 

g{y,r) = ^ ^ (15) 

and consider the derivative of pij) and (|12)) with respect to y we have: 
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dtjy) 
dy 



^ dy 



and 



dy 



2y 



2 ^ 



Tiin \ / 



dr 



dr 



'^giy,r) 

fir) 



drr, 



ygiy,r) 



dr 



dr 



drr. 



(16) 
(17) 

(18) 
(19) 



Comparing the integral terms in the two equations above, we see that they are 
itical upto a factor of y. We 
to rewrite the first equation as: 



identical upto a factor of y. We can then use the fact that at r = rmin, y"^ ~ '^^^""'"^ 



dy 



{9{y,r)f f'2y9{y,r) 

2 ar-y ' 



drr, 



lin dy 



(20) 



After checking that the divergent piece of the integral cancels with the divergent 
second term in pOj) (see appendix B) we therefore have that: 



dt{y) d^{y) 

-y- 



dy 



dy 



which can be rewritten as 



(21) 



(22) 



Thus from the plot of the endpoints, by taking the gradient at each point we 
are able to obtain the value of y for that geodesic. Note that at no point in this 
derivation have we used the fact that the spacetime is asymptotically AdS, and it can 
easily be shown that the result holds for any static, spherically symmetric spacetime 
(i.e. any metric of the form of ()38() ).^ Calculations along these lines are also not 
restricted to null geodesies, see |14i for similar derivations involving spacelike ones. 
In the next section we will see how we can further use this to recreate the metric 
function f{r). 

3.2 Reconstructing /(r): Method I 

The method involves iteratively recovering /(r) starting from large radius and work- 
ing down towards r = 0. We are working in an asymptotically AdS spacetime. This 



■^It is also not necessary for the boundary on which the geodesies begin and end to be at infinity; 
the result is valid for any spherical boundary within the bulk. 
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means that as r ^ oo, /(r) — > r"^ /B?' + 1; thus we can say that for r > Vn for some 
r„ which can be arbitrarily large, /(r) ~ r"^ /E? + 1. We set the AdS radius (ii) to 
one, and by considering large enough y = yn such that the minimum radius, rmm; 
corresponds to r„, (fTTl) becomes (setting the initial time, Iq, to be zero): 



tn= r At =2£^^dr (23) 



2 / dr = TT (24) 

(r2 + l)Jl-y2r£+l 



as we would expect, as geodesies remaining far from the centre of the space would 
not "see" the modification and thus behave as in pure AdS. In this case we can solve 
the equation for the minimum radius: 

to get r„ in terms of y„, which is determined from the geodesic endpoints. Thus we 
can determine r„ (and hence /(r„) = + 1). 

If we now consider a geodesic with slightly lower ratio of angular momentum to 
energy, say with yn-i = yn — where e > 0, we have the equation: 

which can be split up as follows (noting that r„,_i < r„): 

tn-i = 2 %^ dr + 2 dr (27) 



f{r) Jr„ f{r 

The second integral can be evaluated by setting /(r) = r^ + 1 as in ((231), and we 
have 

2 f ^^^^ dr = . - 2 arctan (/(T^^CK^^C ) (28) 

We can approximate the first integral by taking a Laurent expansion about the 
point r„_i. This gives, to lowest order: 



rn-l 



9{yn~i,r) 
fir) 



dr 



rn~l) 



-1/2 



f'{rn-i) 

/K-i) 



dr 



2^^ 



rn-l 



(29) 

(30) 
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where we have used that = r^_i//(?'n,-i) and integrated. This can be simphfied 
further by writing^: 

/ [rn-i) ~ (31) 

Thus we have: 

tn-i^ ^^J:~"''-t , +vr-2arctan(^(l-e,)r^.-,L,) (32) 



r-n-l f(rn-l) 



which can be used in conjunction with y'^-i = r'^-i/ f{rn-i) to calculate r„„i and 
/(r„_i), as we know and yn-i from the geodesic endpoints, and r„ and /(r„) 
have already been calculated. 

For general we split the integral into several pieces; the two "end" integrals, 
which we evaluate as in (|28j) and (|29|) , and a series of integrals in the middle which 
can be evaluated using the trapezium rule or similar (see Fig0]). The formula for 
general tn-i is then given by (for i >2): 

=2 r ^%pi^dr (33) 

~ -^n—i + Bn—i + Cn—i (34) 



where 



fir / 2'-n-i + l _ /(T-n-i + l) _ 1 

«»-£<----4%e^-%er) (3a, 



Cn-i = vr - 2 arctan (^^(l - y2_.) r2 -y2_,.^ 



(37) 



For a spacetime of the form of (|lfl|) with monotonic effective potential, we can 
use the above to recover the function /(r) down to r = 0. See section 5 for examples. 



3.3 More general spacetimes 

So far we have restricted ourselves to considering metrics of the form of (^, however, 
we can look to extend this further by considering the most general static, spherically 
symmetric spacetimes, given by metrics of the form: 

ds2 ^ -f{r)dt^ + h{r)dr'^ + r^d^l (38) 
■*Note that this hnear approximation of the gradient is only fine whilst r„ — r„_i is kept small. 
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2 g [yn-i> r] 

f [r] 




Figure 4: A plot showing how the integrand of (|3.S|) . is spht up into two end curves 
and a number of trapeziums in order for the integral to be well approximated. The 
actual curve is shown in black, with the approximations in red. 

where the dr^ coefficient is allowed to be different to the dt'^ coefficient. Equations 
© and 0) for the geodesies remain unchanged, and the third is now given by: 

r^.(^ + ^ ^ 

\h(r) f{r)h{r) h{r)r'^ 

where k = for the null case. We can combine the equations together as before to 
give the integral equations: 




The question now is whether we can use the two equations above along with 
the expression = r^/ /(ro) relating the parameter y to the minimum radius to 
recover both /(r) and h{r) and thus reconstruct the bulk metric from the endpoints. 
This initially appears possible as we haven't explicitly used our expression for (p, 
(|41|) . in our method of the previous section. However, our observation that the 
gradient of the endpoints is given by the angular momentum H22|) suggests that we 
don't have three linearly independent relations. Attempting to continue as before, 
by approximating the spacetime as pure AdS at some large radius and working 
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inwards, the degeneracy soon becomes apparent (see appendix C), and we are unable 
to recover bulk metric information in these more general spacetimes (unless we a 
priori know the relationship between /i(r) and f{r)). 



3.4 Reconstructing /(r): Method II 

There is an alternative method for reconstructing /(r) in a spacetime of the form 
of which makes the degeneracy of the expressions for t and (p more apparent. If 
we start from expression ()21() and integrate over y we get: 

'd,' = /.|d. (42) 
which can be integrated by parts: 

t{y)=y<P{y)- [ <pdy (43) 



and then rewritten using the expression for (p from ()12|) : 

t{y) = yHy)- [ [ , drdy (44) 

We can now reverse the order of integration, and as the function /(r) has no 
dependence on y, integrate over y: 



fir) 



t{y)=y4>{y)+ I -y^x/l-y^^dr (45) 



r 



Taking the initial conditions to be (0o,to) = (0,0), we can say that for any 
endpoint {(j)j,tj) on the boundary. 



dt 



tj) ^ Jr. f{r)9{yj,r. 



■ dr (46) 



where we have renamed the minimum r as rj and used our definition for g{y, r) from 
before. This, coupled with the equation for the minimum r. 



allows the metric function /(r) to be fully reconstructed from the plot of the end- 
points, using a method of approximating the integral similar to that in section 13.21 
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4 Analysing the two methods 



Whilst the previous methods both in principle allow the bulk metric to be fully 
recovered (from /(r„) down to /(O), with r„ arbitrarily large), there are a number 
of factors to take into account when performing the reconstruction in practice, in 
order to obtain a balance between accuracy and computational effort. 

Consider a setup where we divide the range of y into + 1 equal segments, with 
2/0 = and un+i = 1, so that we have + 1 steps in the process of recovering /(r) 
(Note that the first step in the iteration uses yjsr rather than yAr+i, as we want the 
first minimum radius to be finite). 

There are a number of issues with recovering /(r) in this scenario. We firstly 
observe that splitting y up into equal segments means that the minimum radii will 
not be equally separated; there will be a greater number of low radius points (due 
to the shape of T4jj). This initially appears helpful, as we expect most of the bulk 
metric's deviation from pure AdS to be localised close to the centre, and so a greater 
number of data points in this region should improve our estimate of /(r). 

However, as both methods for recovering /(r) involve approximating integrals 
over r, accuracy is only maintained if the step size is kept relatively small for all r; 
whilst starting with y^r as close to one as possible allows /(r) to be recovered from 
as large a radius as possible (as higher y^ means higher r^r), this also means the 
first few steps in r will be unacceptably large. This can be countered by ensuring 
the number of steps is appropriately large, however, this will then result in a longer 
computational time. For example, in the asymptotically AdS spacetime with /(r) 
given by 

if use a step size in y of 0.0005 and take the initial y^ to be 0.9995, the corresponding 
minimum radius r^v ~ 31.5. The minimum radius corresponding to ^at-i = 0.999 is 
only Tjv-i ~ 22.3 however, and such a large jump in the radius causes problems with 
using our approximation of the form of (|32|) . Whilst the original Laurent expansion 
of H29() is still fairly accurate over this distance, our approximation of the gradient 
by 

/'(r„_.) = -_ f^--^ (49) 

^n—i+l 1"n—i 

is not. For the example above, the actual value of f'{r) at r = 22.3 is 44.5, whereas 
the value given by (|49j) is 53.8; a large discrepancy. 

We can still split the yn linearly (in order to keep the majority of the points at low 
radii), however, we choose the initial yjv to be be slightly lower. This enables our first 
step in the radius to be kept small. ^ In our example of H48|) . if we choose y^ = 0.9985 

^This can be checked in practice by calculating rjv and tn-i for a given choice of initial j/jv; if 
the difference rjv — rjv-i is too great, either lower j/jv or increase TV until it's acceptable. 
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Initial y Step size a (4) /3 (1) 7 (8) 



0.9995 0.0005 0.778 0.971 0.971 
0.9985 0.0005 1.39 1.52 1.52 



0.9995 0.00005 2.37 2.24 2.24 
0.9985 0.00005 3.42 1.13 6.18 



Table 1: Comparing the accuracy of using two different values for ^at to recover the 
bulk metric information. The actual values for a, /? and 7, which correspond to 
the three numerical factors in (|48j) (and are properly defined in (|53|) ). are given in 
brackets. More details on how these estimates were generated is given in section 5. 

and keep the step size at 0.0005, we find that r^v — ?'Ar-i ~ 2.5 which whilst still 
quite large, is much more acceptable than the value of 9.2 we had previously. Table 
n shows how this discrepancy in the gradient affects the extraction of /(r) by the 
method of section for the example of H48|l : it is worth noting at this stage that 
even with the lower choice of starting y the recovered estimate for /(r) is still quite 
poor with a step size of 0.0005. This is important because although we could lower 
j/AT even further, we must keep our initial radius reasonably large; if we choose 
to be too low, we risk having being too small, such that our assumption that the 
spacetime there is approximately that of pure AdS is no longer valid. ^ We can keep 
our choice of y^ high by reducing the step size (as in the lower half of Table but 
at the expense of additional computational time (For more details on the numerical 
extraction see the examples in the next section). 

There is, however, another problem with numerically reconstructing /(r) using 
the approximation of the gradient by (|49j) . If we examine the form of the term 
containing f'{r), namely the An-i in H33|). whose equivalent in method II is: 



we note that the f'{r) appears with a negative sign inside the square root. As our 
approximation of f'{r) by (|49|) is in general an overestimation of the gradient, as for 
large r, /(r) oc + 1, when solving for r„_j and /(r„_j) the term under the square 
root can become negative, leading to the estimates being imaginary. This can occur 
even if the overestimation is small. 

There are a number of possible ways in which to resolve this problem. We can 
firstly consider a different approximation of the gradient to the one given in (|49j) : 

^This poses an interesting question as to how one determines what a "reasonable" initial radius 
is for an unknown spacetime. If our rjv was too low, such that the metric was not approximately 
pure AdS at this point, how would this be apparent from the extracted estimate of f{r)? A simple 
solution is to check the behaviour of the estimate at r close to rjv; if the function does not continue 
as r'^ + 1 (i.e. pure AdS) down to say r = rjv/2, then the chosen initial radius was too small. It 
is also worth noting that in a physical situation, one would naturally expect ordinary matter to 
remain within a radius of order of R from the centre, due to the confining AdS potential. 




(50) 
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the over-estimation can be avoided by using an expression such as 



J [rn-i) = ^ ( + ) (51) 

which takes the average of the two nearest hnear tangents to the curve. Inserting 
this expression into either of our iterative methods for recovering /(r), however, 
immediately raises a new problem from a computational point of view; we now 
cannot recover a value for r„_j and /(r„_j) at each step, as the n — i term now also 
depends on the subsequent term, n — i — 1, as well as the previous. We now have 
to determine all the equations for the various steps and solve them together at the 
end to recover f{r)J This is a much more complicated operation than numerically 
solving at each step, and leads to a considerable increase in the computational effort 
required to reconstruct the metric function. 

Another approach is to take a higher order series expansion around the minimum 
radius in our approximation of the integral, however we then need an expression for 
the second (or higher) derivative of f{r), which leads to an escalation of the problems 
mentioned above; either the overestimation caused by only using approximations 
involving larger r terms upsets the recovery of r„_j and f{rn-i), or we need to wait 
and solve all the equations at the end of all the steps, again greatly increasing the 
computational effort required. 

The gradient of /(r) appears in the integral of because of the Laurent 
expansion we have used in its approximation; in method II, however, we can use a 
different approximation which avoids this. 

FigEl shows a plot of a typical curve we are trying to approximate. The crucial 
difference between this type of curve and that of Fig^ is the behaviour at the 
minimum radius, which is due to the positioning of the g{y, r) term in equations 
(|33|) and ()46() . That g{y,r) appears in the denominator of H46|) rather than the 
numerator results in the curve of FigO having an infinite gradient at the minimum 
radius. We can then describe the leftmost part of the curve as following a vertical 
parabolic path. Thus instead of using a Laurent approximation at this point, we 
can make a geometric one and say that: 




4 (r„_i+i - rn-i) L y 



,2 



An-i = - — J- ^ — xr--— — ^2) 

3 /(r„_i+i) y 



where we have used that the area under a parabola is 2/3 the width multiplied by 
the height. The gradient of the curve will always be infinite at the minimum radius, 
and as the other terms {Bn-i and Cn-i) do not depend on f'{r), we avoid the need 
for an approximation to the gradient of the form of l|49jl. Importantly, we also avoid 
any increase in the computational effort required to recover f{r). It is also worth 
noting that there is no similar geometric argument for such an approximation to be 
applied in the original method of section 1221 For similar step size this modification 

^Note that the final gradient f'{ro) will have to be approximated using an expression of the 
form of lUni 
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Figure 5: A plot showing an example curve (the integrand of (|4()j) ) to be approxi- 
mated in method II. The leftmost part of the curve is initially a vertical parabola. 

leads to a considerably more accurate estimation of /(r) from the endpoints, as we 
shall see in the next section, where we use both methods to reconstruct /(r) in two 
different cases. 

5 Examples 

Consider a small deviation from pure AdS, such as the spacetime with metric func- 
tion /(r) as in (|1H|) . Choosing our initial value of y to be = 0.9985 (for the 
reasons of the previous section), we investigate the ability of the two methods to 
extract the bulk information. FigEl shows our approximations of /(r) generated 
using method I for a range of different step sizes. 

For this first example, the top set of curves in FigEl suggest that we are easily 
able to extract the bulk metric using any of the step sizes, as they all appear to lie 
very close to the actual function /(r). However, this is somewhat misleading, as we 
shall see. In order to more closely examine the accuracy of the recovered bulk data, 
we use a non-linear fit of the form: 



in Table [21 As we mentioned above, looking at the top plot of FiglHl all of the 
approximations seem very close to the original function, however. Table |21 shows 
this not to be the case. Until we go to a very small step size, we are unable to 
accurately extract what one might consider the important metric information; if the 




(53) 



to obtain values for a, (5 and 7 for each estimate, and the results are presented 
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Initial y 


Step size 


a (4) 


P (1) 


7 (8) 


0.9985 


0.002 


0.371 


0.530 


0.530 


U.99oO 


O.UOi 


0. /9 ( 


0.995 


0.995 


0.9985 


0.0005 


1.39 


1.52 


1.52 


0.9985 


0.0002 


2.18 


2.11 


2.11 


0.9985 


0.0001 


2.74 


1.47 


3.92 


0.9985 


0.00005 


3.42 


1.13 


6.18 


0.9985 


0.00002 


3.80 


1.11 


7.20 



Table 2: Best fit values (to 3 s.f.) for a, /3 and 7 for data generated using method I, 
with the actual values indicated in brackets. Only in the lower half of the table do 
the estimates for the three unknowns really start to converge to the correct values. 



Initial y 


Step size 


a (4) 




7(8) 


0.9985 


0.002 


1.61 


2.01 


2.01 


0.9985 


0.001 


3.34 


1.16 


6.01 


0.9985 


0.0005 


3.92 


1.01 


7.77 


0.9985 


0.0002 


3.99 


1.00 


7.97 


0.9985 


0.0001 


4.00 


1.00 


8.00 



Table 3: Best fit values (to 3 s.f.) for a, (3 and 7 for data generated using method 
II, with the actual values indicated in brackets. Using the alternative method, we 
are able to pick out rough values for the three unknowns as early as with a step size 
of 0.001, and by a step size of 0.0001 the estimates have converged to the correct 
values (to 3 s.f.). 

modification of the spacetime caused by the extra term in (|48|) was to correspond to 
some physical phenomenon, we might expect its "mass" and "extent" to be related 
to the quantities a, and (3 and 7 respectively. Whilst the best fits we obtain are 
converging to the correct values as we take the step size smaller, we are already 
calculating a significant number of terms to generate these estimates. 

If we compare this to the fits we obtain using the data calculated via method II, 
with the modification to the A^-i terms given in (|52|1 , we see a notable difference in 
the results. The estimates again all look good when plotted with the actual function 
/(r), as we see in Fig|3 however, as the bottom plot shows, this accuracy is now 
maintained at smaller scales. Table |31 shows the new values for a, (3 and 7 for each 
step size, obtained using the same non-linear fit, (|53|) . as before. 

Comparing these estimates with those in Table |21 we see that the recovered 
values for a, (3 and 7 are far more accurate using this alternative method, for each 
choice of step size. We now have a very good fit to the actual values of 4,1 and 8 
in obtained using a relatively low number of steps, and thus relatively quickly. 
Indeed, we see that to obtain the same degree of accuracy using method I we need 
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Figure 6: Plots showing the various estimates for /(r) given in Table|2l(i.e. produced 
using method I), compared to the actual metric function /(r) from ()48() . Whilst all 
the estimates seem good fits over a large radius (top figure), closer consideration of 
the curves highlights their differences (bottom figure), which lead to the inaccurate 
values calculated for a, (3 and 7 from the larger step size estimates. 
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Figure 7: Plots showing the various estimates for /(r) given in Table |21 compared 
to the actual metric function /(r) from (|48j) . Using method II for producing the 
estimates has resulted in much closer fits, both over a large radius (top figure), and 
at smaller scales (bottom figure). This is refiected in the highly accurate values of 
a, P and 7 determined from the all but the largest step size estimates. 
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Initial y 


Step size 


a (4) 


/? (1) 


7 (8) 


X (3) 


V (2) 


A (1) 


0.9985 


0.002 


-1380 


133 


133 


2.60 


2.20 


1.20 


U.99oO 


U.UUi 


U.Uoz4 


U.4ZD 


U.42D 


2.69 


z.iz 


i.U4 


0.9985 


0.0005 


0.790 


0.880 


0.880 


2.96 


2.03 


0.981 


0.9985 


0.0002 


1.69 


1.25 


1.25 


3.06 


1.95 


0.940 


0.9985 


0.0001 


2.23 


1.62 


1.62 


3.05 


1.94 


0.946 


0.9985 


0.00005 


2.87 


2.39 


2.39 


3.00 


1.97 


0.985 


0.9985 


0.00002 


3.68 


1.31 


6.50 


3.09 


2.00 


1.04 



Table 4: Best fit values (to 3 s.f.) for a, 13, 7, Xi V ^ data generated using 
method I, with the actual values indicated in brackets. In this more complicated 
modification to AdS, the presence of the sine term masks the finer structure and 
prevents the original method from converging towards the correct values of (3 and 7 
until our smallest choice of step size. 



to use step sizes which are more than twenty five times smaller, which (assuming 
a standard computational time per step) would therefore take at least twenty five 
times longer to compute. Thus using the modified second method offers a significant 
improvement over the original. 

Can we use the methods to recover a more complicated /(r)? If we consider a 
spacetime similar to that in (|48j) . but which is further modified at low r by an extra 
term: 

^/ N 2 3rsin(2r) , ^, 

we can use both methods to generate estimates for f{r), and using a non-linear fit 
of the form of: 

nr)-l+r (^2 + ^)(^2+^) + ^4 + ;, 1^5) 

we can judge the accuracy of the estimate against (|54|) . Using the same starting 
value for y and the same step sizes as before, we obtain the results shown in tables 
m and 13 and figures |H1 and El 

Once again we see that the second method produces a much better estimate of 
/(r) than the first; the difference in the accuracy is again significant. What is also 
noticeable in this slightly more complicated example is that the presence of a more 
prominent feature in the metric function (in this case the small oscillation from the 
sine term) can obscure the finer details of the spacetime. This is most obvious in the 
results from using method I (see Table |1J where we can clearly see that whilst the 
numerical factors corresponding to the sine term (x, rj and A) are easily determined 
even at a large step size, it then requires a smaller step size to pick out the a, f3 and 
7 values than it did in the first example. 
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Figure 8: Plots showing the various estimates for /(r) given in Table ^ (i.e. pro- 
duced using method I), compared to the actual metric function /(r) from (|54jl . The 
convergence of the estimates to the actual curve can be seen in both plots, and 
whilst the values obtained for Xi f] and A are reasonably good for all step sizes, only 
upon using the smallest step size do we start to pick out the a, f3 and 7 values. This 
is why, in the lower figure, the red curve representing the final estimate appears to 
be a worse fit to the actual curve than some of the others; the range of r shown 
highlights the accuracy of the Xj V and A values rather than the a, /3 and 7 ones. 
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Initial y 


Step size 


a (4) 




7 (8) 


X(3) 


77 (2) 


A(l) 


0.9985 


0.002 


1.17 


1.46 


1.46 


3.47 


1.94 


1.12 


0.9985 


0.001 


2.90 


2.52 


2.52 


3.03 


1.97 


0.998 


0.9985 


0.0005 


3.85 


1.06 


7.38 


3.00 


2.00 


1.00 


0.9985 


0.0002 


3.98 


1.01 


7.93 


3.00 


2.00 


1.00 


0.9985 


0.0001 


4.00 


1.00 


7.99 


3.00 


2.00 


1.00 


0.9985 


0.00005 


4.00 


1.00 


8.00 


3.00 


2.00 


1.00 



Table 5: Best fit values (to 3 s.f.) for a, /3, 7, x-, V aiid A for data generated 
using method II, with the actual values indicated in brackets. Once again, the 
alternative method proves much more adept at accurately estimating the unknowns, 
and converges to the correct values almost immediately. 




Step size 

0.002 

— 0.001 

0.0005 

0.0002 

0.0001 

0.00005 

Actual 



0.7 0.8 



Figure 9: Plots showing the various estimates for /(r) given in Table|Sl(i.e. produced 
using method II), compared to the actual metric function /(r) from (jHU). Only the 
estimates generated using the largest two step sizes are visibly distinguishable from 
the actual curve, as the the values of q, /?, 7, X) V and A are all quickly picked out 
using this alternative method. 
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Whilst the second example has shown that the methods can be used for /(r) 
which have a significant deviation from pure AdS near the centre, there remain 
problems in applying the methods in certain circumstances, which are examined in 
the next section. 

6 Limitations 

There are limitations in recovering /(r) all the way down to r = in the cases where 
the metric is non-singular but does not give a monotonic effective potential for the 
null geodesies, and also when there is a curvature singularity at r = 0. 

6.1 Metrics with a non-monotonic effective potential 

There is a problem with the method for reconstructing /(r) if the metric is such 
that for some critical value of the angular momentum, the null geodesies can go 
into an unstable orbit (assuming fixed energy). This is indicated on the endpoint 
plot by the <j) coordinate heading off to infinity, and the presence of a jump in the 
time taken by the geodesies with L < Lcrit compared with those with L > Lcrit, see 
Fig llUl for an example. 

The reason for these unstable orbits is because the effective potential for the 
metric isn't monotonic; the orbits occur when the greatest radius (r^ say) for 
which dV/dr = corresponds with that for which 1/ = 0, see Fig llll This non- 
monotonicity also explains the time delay for those geodesies with angular momen- 
tum lower than Lcrit, indicated by the gap between the asymptotes in FiglTIll those 
with slightly lower L have an effective potential of the form shown as the dashed 
blue curve in Fig llll and travel down to minimum radius whereas those with 
L > Lcrit have minimum radius greater than or equal to r^, indicated by the dashed 
red curve. The extra (non-negligible) distance Vp — fp—i travelled by the geodesies 
results in the time delay seen in the endpoints. 

We thus have the situation where we have no geodesies with minimum radius 
rq, where rp_i < < r^, and this is where our iterative methods break down. 
Although we can begin by using one of the methods from earlier to recover /(r) 
down to r = Tp, which corresponds to using the endpoints in the "lower branch" of 
Fig llUl we cannot continue past this point. Using method I, say, the final term of 
the iteration will be: 

r°° 1 

tp = 2 ^=^= dr^Ap + Bp + Cp (56) 

with Ap, Bp and Cp defined as in section 13.21 As the next term would involve an 
integration over the range rp_i to rp, which cannot necessarily be taken to be small, 
it cannot be well approximated by the methods we have used previously. 

The gap between the upper and lower branches of Fig llUl tends to a constant 
time delay, t^, as — > oo, and it is this which corresponds to the integral we are 
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t 

10 r 




3.5 4 4.5 5 5.5 6 

Figure 10: Plot of the endpoints of null geodesies which all begin from the same 
point on the boundary and pass through an AdS-like space time with non-monotonic 
Veff. The red lines indicate the limits to which the black curves tend as L ^ Lcru 




Figure 11: Plot of the effective potential for three null geodesies: one with the 
critical angular momentum, Lcrit (the solid black curve), one with L < Lcrit (the 
lowest curve, dashed blue), and one with L > Lcrit (top curve, dashed red). 
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unable to well approximate: 



rrp I 

td = 2 dr (57) 

where we note that in the limit ^ oo we have that i/p-i = Up. 

We are left with the situation where we do not have enough information to carry 
on recovering /(r) past the rightmost maximum in the effective potential; there is 
no approximation we could attempt to make without already knowing more about 
the form of V^//. Although we have the extra endpoint data (i.e. the top branch of 
EH, we are unable to use it to extract f(r). 

Having discovered this apparent limitation in the extraction of the bulk metric, 
we should ask ourselves how often we would expect to find metrics which give non- 
monotonic effective potentials for the null geodesies. In other words, how often do 
we find spacetimes which allow stable orbits of light rays? Whilst it is very easy to 
construct metrics by hand which allow this, in physical systems circular orbits of 
light rays are a rarity. Indeed, recent work by Hubeny et al., [22], indicates that for 
metrics corresponding to a gas of radiation in AdS, it is impossible for circular null 
geodesies to occur. 

Admittedly our restrictive case of metrics of the form of are unrealistic, as 
any metric of this form necessarily has opposite sign pressure and energy density 
components of the stress tensor, and as our methods currently do not work in the 
more general metrics, ()c{8() . this is merely an observation. 

One scenario in which null geodesies can go into circular orbits is around black 
holes, and we now conclude the current section by considering such an example. 



6.2 Metrics containing a singularity 

We ask what information can be retrieved if the metric contains a singularity at 
its centre, for example Schwarzschild-AdS. We again cannot recover the entirety 
of /(r), as the geodesies which pass behind the horizon will not return out to the 
boundary. Thus whilst in the non-monotonic case we had extra information (the top 
branch of Fig llflj) available but no way to get past the jump in the minimum radii 
of the geodesies, here we only have a reduced spectrum of endpoints from which to 
recover /(r). 

The endpoints we do obtain are those with angular momentum above some 
critical value, and so we can use the above iterative method to recover /(r) down to a 
certain critical radius (which will be greater than the horizon radius). However, any 
null geodesies which pass behind the horizon will always terminate at the singularity, 
and so we are unable to obtain any further information using our iterative method.^ 

*As mentioned in the introduction, there are numerous papers devoted to the extraction of 
behind the horizon information, and we do not concern ourselves with this here. 
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For an example we consider five dimensional Schwarzschild-AdS, in which the 
metric function /(r) is given by: 

where R is again the AdS radius, and rh is the horizon radius. We set i? = 1 as 
usual, and note that the effective potential for null geodesies is given by: 

For a null geodesic to avoid hitting the singularity at r = 0, we require 14// = 
for some r > r^. The lowest L for which this will happen will be when the gradient 
of Veff is also zero. Differentiating (|59|) with respect to r, and setting equal to zero 
gives 




which is equal to zero when r = ^/2^Jrf^ + r^. Thus for the case where = 1 say, 
the minimum r which can be probed by the null geodesies is r = 2, and hence we 
can obtain information about the metric function /(r) from r = oo down to r = 2. 



7 Discussion 



We have seen that the endpoints of null geodesies can be used to recover information 
about the bulk in asymptotically AdS spacetimes. By varying the ratio of angular 
momentum to energy of the geodesies we are able to systematically probe to different 
radii, and this, combined with the appropriate AdS boundary conditions, allows us to 
reconstruct the metric. What is perhaps quite surprising is that this local, geometric 
information about the spacetime can be completely recovered in certain scenarios; 
in a metric of the form of , as long as the effective potential for the null geodesies 
is monotonic, we are able to completely recover the function f{r). Whilst this is 
of course due to the fact that we have severely restricted our spacetime in such a 
scenario, to ensure that the amount of information we require is contained within 
the bi-local spectrum of endpoints, it is still impressive that the recovery process 
itself is so simple. Its simplicity lies in the remarkable nature of our expression for 
the gradient of the endpoints of null geodesies 1)22(1 : to have that the relative angular 
momentum of the geodesic corresponds to the endpoint in such an accessible fashion 
is astounding, and forms the basis of our ability to reconstruct the bulk. 

Focusing on the cases with a monotonic Veff, we developed two different meth- 
ods for carrying out the recovery of the bulk information, both of which involved 
iteratively reconstructing /(r) from large r down to the centre. In the process of 
actually using these methods to numerically extract the bulk information, we made 
several refinements to increase their efficiency, the most important of which involved 
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our approximation of the gradient of /(r). In both our original methods, we were 
consistently overestimating /'(r), a problem that was significantly affecting their 
accuracy, unless we kept the iterative step size very small. Whilst we suggested a 
number of ways of resolving this problem, all made the process of recovering /(r) 
considerably more complex to implement, thus requiring a great deal of extra com- 
putational effort. Given that one of the main aims in this work was to produce a 
practical method for the recovery of the metric information, these resolutions were 
not acceptable. Rather than investigating further methods of approximating f'{r), 
we looked at the question of why we needed an approximation to the gradient of 
/(r) at all, as f'{r) does not appear anywhere in our original geodesic equations. 

The solution lay in finding an alternative to the Laurent series we were using to 
approximate the section of the geodesic path close to its minimum radius. Whilst 
this had been the natural approximation to use in our first method, there was a 
simpler one we could apply in the second, which stemmed from a crucial difference in 
the type of curve we were trying to approximate. In method I, the curve in question 
is infinite at the minimum radius, because the g{y,r) term (see (|15|)) appears in the 
numerator of ()33() . In method II, on the other hand, the g{y, r) term appears in the 
denominator of the relevant equation, H46() . which results in the curve being zero with 
an infinite gradient. This allowed us to approximate the curve geometrically by a 
vertical parabola, and thus the integral by the corresponding parabolic area formula. 
This modification enabled us to estimate /(r) far more efficiently, as it removed any 
reference to f'{r). It also greatly improved the ability of the iterative process to 
extract what one might consider the important metric information: the numeric 
values relating to the "mass" and "extent" of whatever physical phenomenon was 
causing the deviation from pure AdS. 

We concluded the project by considering two alternative scenarios: where the 
effective potential was non-monotonic, and where the metric contained a singularity 
at the centre. In each case, we were only able to recover /(r) down to a certain 
critical radius, corresponding to a local maximum oi V^ff (see fig Jllj) . We were 
unable to continue the process of recovering /(r) any further as, in the non-singular 
case, there is a non-negligible "jump" in the minimum radii of the geodesies, due 
to the shape of the effective potential, which we cannot well approximate. Thus 
although we have additional endpoint data, we are unable to use it to extract more 
of the bulk information; we do, however, note that in physical situations, which is 
where we would eventually hope to apply these techniques, we do not often encounter 
such metrics (i.e. spacetimes in which light rays can enter into circular orbits). For 
those metrics containing a singularity, we do not have any further endpoints with 
which to continue recovering /(r), as any null geodesic travelling to lower radii will 
pass behind the horizon and terminate at the singularity. 

A possible extension of this work is in considering less restrictive forms of the 
metric, such as those where we have relaxed spherical symmetry, or where the met- 
ric can evolve with time. Whilst our previous attempt (in section 3.3) at modifying 
the iterative methods to work in a more general spacetime failed due to the extra 
unknowns we introduced, this does not mean every modification would be so prob- 
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lematic. One can envisage, for example, a situation in which the metric function 
fluctuates over time, in such a way that the effective potential for the null geodesies 
remains monotonic throughout. In this scenario, whilst the process might be more 
complicated, we should still be able to extract the entire bulk metric, including the 
time dependence, from the complete set of null geodesic endpoints. 

Another natural extension of this project would be to consider the endpoints of 
spacelike geodesies, and examine whether they can be used to recover more infor- 
mation about the bulk in the various scenarios. The main issue with producing a 
method to iteratively recover the metric information as we have done here with the 
null geodesies is that our equation (|21j) involving the gradients of the endpoints is 
significantly more complicated. For spacelike geodesies we have two independent 
parameters (L and E) to consider; simply using their ratio as we did in the null case 
is no longer sufficient. Other works have examined this in more detail, e.g. and 
they give encouragement that progress can be made in this direction. The challenge 
then would be whether a computationally practical method can also be found. 
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Appendix A: Null geodesic paths in AdS space 

For the null geodesic paths we have x'^ = 0, and we can solve the following three 
equations numerically to provide the plots of the geodesies: 

L = r'^ip (62) 



2\ ^ ^2 



1 + -^ *' + T^ + '-V=0 (63) 



It is also useful to calculate the geodesies analytically. This can be done most 
easily by rewriting the original equation for the metric Q in a new form. Rescaling 
the radial and time coordinates r — > i?tan(r) and t ^ Rt gives: 

= _i?2 (gec2(r)dt2 + sec^ (r)dr2 + tan2(r)#2) (54) 

where we have suppressed two of the angular coordinates as before. This leads to 
the modified constraint equations: 

E = R^sec^{r)i (65) 
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(66) 



sec\r)t^ + sec^ (r)r^ + tan^(r)(/.^ = (67) 



Combining the three equations, by substituting and into (|H7j) to ehm- 
inate t and (/) we have: 

+ sec^ (r)r2 ^ = o (68) 



i?^sec^(r) i?2tan^(r 
Rearranging, and again using H65|) to ehminate the dependence on A gives: 



^2 . dr dr E 



(i?2sec2(r))^ i?4tan2(r)sec2(r) dA dtR^sec^{r) 



and so 



Using the substitution x = ^"'^(0 ^j^jg equation becomes: 

dx 



which imphes x = — sin(t + t^) and thus that 



We can also calculate the dependence of </> on t: 

L ■ dS dt dd) E 



i?2tan2(r) ^ dA dtR^sec^{r) 

thus: 



(69) 



dt~y E^sin^r) ^ ^ 



dt 



^2 

cos(r) = - ^/ 1 - — sin(t + to) (72) 



(73) 



(74) 



dt Esin^ir) 
which can be combined with 1)72^ to give: 

tan{(l) + (po) = ^ tan{t + to) (75) 
E 

Note that the dependence of the geodesies on the AdS radius of curvature has 
been scaled out in ((7^ and (|7S)) by the choice of coordinates. 
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Appendix B: Cancelling the divergent term from the 
Leibniz rule 

To make sure our expression for the derivative of t with respect to y is finite, we 
have to combine the two terms from: 



dtjy) ^^d_ 
dy dy 

= 2y 



^ 9{y,r) 
fir) 
{9{y,r)f 



dr 
dr 



'^giy,r) 

fir) 



drn 



1171 

dy 



(76) 
(77) 



to ehminate the divergence piece. We need an expression for drmin/dy, obtained 
using our expression for y in terms of the minimum radius, (|47)1 . 



dr„ 



dy 



dy_ 

dr 



f'ir„ 



' mm J y mm 



2 firmin) 



(78) 



which we then use to rewrite the second term of (|77|) as 3,n intcgrcil from r^^^ to 



oo: 



( My, 



V fir) 



drr, 



dy 



if'ir)-rf"{r))^W)9iy,r) 



ifir)-y'ir)y 

fir)g{y,r) 



dr 



+ 



ifir) - i/'(r)) /7(r 

2y (9(y,0)' 



dr 



dr 



We see that the third term above cancels precisely with the divergent integral 
in H77|) . and after rearranging we are left with 



dtjy) 
dy 



r{f'{r))^-2rf"{r)f{r)g{y,r) 



2(/(r)-§/'(r))^/7H 



dr 



(79) 



which is finite. 



Appendix C: Attempting to recover the bulk metric in 
the more general spacetime 

In the more general case of a spacetime described by equation (|38|) . we have two 
unknown functions /(r) and /i(r) to recover, along with r itself, and hence we need 
three linearly independent equations. We appear to have them in equations ()4U() . 
(jH} and (EH. 
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Following the method of section l3?2l we take r large enough so that the spacetime 
is approximately pure AdS, and solve 



= ttS = ^ (80) 

to give r„, and thus also /(r„) = l//i(r„) = + 1. 

We now consider a geodesic with slightly lower y and follow the method as before, 
except that we now use the expressions for t and (p given in (|4()j) and (|41j) . Beginning 
with the expression for t, we split it up as in ()27() and evaluate the second integral 
as in ((2HI)- This leaves 

t„_i « 2 / " ^ dr+7r-2arctan f ^ (l - r2 - (gi) 

ir„_i /(r) / 1 _ V / 

where we can approximate the first integral by its lowest order Laurent expansion 
around the minimum radius. This gives 



~ . + vr - 2 arctan a/ " yn-i) - Vn-i (82) 



Similarly for the expression for (;/) we obtain 



« J^^^^^^^^ + 2 arctan ^"'^ I (83) 



which initially appears to be linearly independent of (|82j) . as although the first term 
in (|83() is simply the first term in (|82() divided by yn-ij the two second terms appear 
different. Even rewriting the second and third terms of (|82|) using 



^ - arctan (j (l - ) ^2 _ y2_ J ^ g^^.^^.^^ | i (84) 

does not linearly relate to the second term of (|83|) . as there is the extra Un-i factor 
inside the arctan. 

The degeneracy occurs because we are attempting to use these two equations, 
along with yl_^ = rl_^/ f{rn-i), to recover r„_i, /(r„_i) and /i(r„_i); we already 
have values for the other variables. Thus the only term which is of importance 
in recovering the unknowns is the first, and the apparent independence of the two 
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equations has arisen from our approximation of the second term. Hence equations 
and (IHSl) effectively both reduce to: 

^ " ( i) _ (.Qfigiant (85) 

We are thus left with only two equations to determine the three unknowns, and so 
the method of extracting the metric functions immediately breaks down.^ 
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